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We propose in this paper a scenario of spontaneous baryogenesis in cosmological models of 
Quintessence by introducing a derivative coupling of the Quintessence scalar Q to the baryon current 
Jg or the current of the baryon number minus lepton number Jg_g- We find that with a dimension- 
5 operator di_iQJg_g suppressed by the Planck mass Mpi or the Grand Unification Scale Mgut, 
baryon number asymmetry ub/s ~ 10“^° can be naturally explained via leptogenesis. We study 
also the isocurvature baryon number fluctuation generated in our model. 
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Evidence is increasing that Universe is spatially flat and accelerating at present time Q . The simplest account of 
the cosmic acceleration seems to be a remnant small cosmological constant A, however many physicists are attracted 
by the idea that a new form of matter is causing the cosmic acceleration. This new form of matter is not clustered 
gravitationally on the scale of the galaxy clusters and has been usually called dark energy or Quintessence |||-^. 

For Quintessence, typically one assumes the existence of a scalar held Q (or multi scalar fields) which is taken to be 
homogeneous in space and via its kinetic and potential energy contributions to the energy-momentum tensor tuned 
to provide an equation of state leading to accelerated expansion, beginning to dominate the matter content of the 
Universe today. Models of Quintessence can have interesting cosmological properties including tracking behavior. 

Generally minimal or non-minimal couplings of Quintessence with gravitation are considered in the literature 
Being a dynamical field. Quintessence is expected to have interactions with the ordinary matters 0, however as 
argued in Refs. for most of cases the couplings are strongly constrained. But there are exceptions. For example, 
Carroll Q has considered an interaction of form with being the electromagnetic field strength tensor 

and studied its implication on the rotation of the plane of polarization of light coming from distant sources. In this 
paper we introduce a type of interaction of Quintessence with the matter, which in terms of an effective lagrangian is 
given by 

^eff = ( 1 ) 

where M is the cut-off scale which for example could be the Planck mass Mpi or the scale of Grand Unification Theory 
Mgut, and c is the coupling constant which characterizes the strength of Quintessence interacting with the ordinary 
matter in the Standard Model of the electroweak theory. Specifically we take in this paper to be the baryon current 
Jg or the current of baryon number minus lepton number Jg_g, and study their implications on baryogenesis. The 
lagrangian in Eq.(|^ involves derivative and obeys the symmetry Q ^ Q constant [^, so the Quintessence potential 
will not be modified by the quantum corrections. 

The mechanism of generati^ the baryon number asymmetry in this paper follows closely the spontaneous baryoge¬ 
nesis |H,|T^. The term in Eq.(P, when Q is non-zero during the evolution of spatial flat Friedmann-Robertson-Walker 
Universe, violates GPT invariance and generates an effective chemical potential fj.^ for baryons, i.e., 


-^dpQJ^ c nj,) , 

Q 

In thermal equilibrium the baryon number asymmetry is given by (when T rrih) 

ObT'^ ,Pb , gbQT'^ 

where gt counts the internal degree of freedom of the baryon. Using the familiar expression for entropy density 

" = , 


( 2 ) 


( 3 ) 


( 4 ) 


we arrive at the final expression for the baryon to entropy ratio 
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Q in Eq.(||) can be obtained by solving the equation of motion of Quintessence given below 

Q + 3HQ + V'{Q) = + SHub) , 


( 5 ) 


( 6 ) 


where H is the Hubble constant and V(Q) is the potential of Quintessence field. For the radiation dominated era the 
Hubble constant is 


H = — = 1 . 665 */^— . 

2t Mpi 


( 7 ) 


The right-handed side of Eq.(||) is about + HQ) and can be neglected unless in the very early universe 

when the temperature T is close to the cut-off scale. 

In this paper we consider a model which has the tracking property, where the potential has a modified exponential 
form 0 , 




V{Q) = fiQ)e 


( 8 ) 


tQ ■ 


In Eq.(^, f{Q) is a function of Q which is required to change more slowly than e "p' in the regime of tracking. 
For the detailed discussions on the specific form of f{Q) and model properties we refer to Ref. [0. Here we simply 
summarize the salient features of this model. For > 3{1 + wb) with wb being the ratio of the pressure to the 
energy density of the background fluid, this model has an attractor solution | 0 - 0 , 


fig = — ^ — = ^2(1 + Wb) 
PQ + Pb 

WQ = Wb ■ 
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During radiation dominated era, pb = and wb = 1/3, so we have 




PQ - 
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Thus we obtain the ratio of baryon number to entropy, which is given by 


ub 5 c 

s “ 7rv^lO(A2-4) 


-1/2 T 
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Taking A^ > 20 


constrained by nucleosynthesis, and 5 * ~ 0 ( 100 ), gb ~ 0 ( 1 ) we have 

s M 
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In the scenario of spontaneous baryogenesis, the baryon asymmetry is generated in thermal equilibrium. This 
requires baryon number violating interactions occur rapidly (T^ > H). However if the B-violating interactions keep 
in equilibrium until Q —> 0, the final baryon asymmetry will be zero. Denoting the epoch when the B-violating 
interactions freeze out by To, i.e. T^{Tb) — H(Tb)), the final baryon number asymmetry obtained is 



0.01 C 

M 


( 13 ) 


For Td around the Grand Unification scale lO^^GeU, Eq.(0) shows that it is quite easy to have ub/s ^ 10“^° 
required by nucleosynthesis and measurement of Cosmic Microwave Background anisotropy, for both M = Mpi and 
M = Mgut- However, if the B-violating interactions of Grand Unified Theory conserve B — L, the asymmetry 
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generated will be erased by the electroweak Sphaleron [^. In this case Sphaleron B + L interactions will make To 
as low as around 100 GeV and ub/s generated will be too small unless c is unreasonably large. 

Now we discuss the mechanism of leptogenesis |l^ in our model by identifying with Jb-l Eq.(|^). Doing the 
calculation with the same procedure as above for we have the final asymmetry of the baryon number minus 

lepton number 


nB-L 
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0.01 c 


M' 


(14) 


The asymmetry ub-l in ([l^ ) will be converted to baryon number asymmetry when electroweak Sphaleron B + L 
interaction is in thermal equilibrium which happens for temperature in the range of lO^GeV ^ lO^^GeV. To in ( p^ ) 
is the temperature below which the B — L violating interactions freeze out. 

In the Standard Model of electroweak theory, B — L symmetry is exactly conserved, however it could be violated 
by the (majorana) neutrino masses via higher dimensional effective operators. For instance, we take a dimension-5 
operator. 


Lf, = -jlLlL(t>4’ + h.c., 


(15) 


where / is a scale of new physics beyond the Standard Model which generates the B — L violations, and (j) are the 

left-handed lepton and Higgs doublets respectively. When the Higgs field gets vacuum expectation value < (j) v, 
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the left-handed neutrino receives a majorana mass rriy ~ 

In the early universe the lepton number violating rate induced by the interaction in (|^) is jl^] 


rp3 

0.04 


(16) 


This rate at To is required to be slower than the Universe expansion rate ~ 1.66gl^^T^/Mpi, which will put a lower 


limit on the scale /, consequently a upper limit on the neutrino mass rrii, |19|. 

To get the upper limit on m^, we try to minimize the value of /. This can be done by minimizing the freezing 
out temperature Tb- From Eq.(14), one can see that if requiring c < 47r so that the effective lagrangian (Eq.(l)) 
is sensible for the perturbative calculation, Tjy must be larger than 10^° GeV for M = Mpi, and lO^GeV for M = 
Mgut ~ lO^^GeV. As a result of it, / should be larger than about 10^^ GeV and the neutrino mass should be smaller 
than ~ 4 eV for M = Mpi. In the case that M = Mgut, f > 4.9 x lO^^GeV and rrii, < O.lkeV. And for both cases, 
Td are much smaller than the cut-off scales. This makes our approximation of neglecting the r.h.s of Eq.(6) valid. 
We have also done a numerical calculation with result which supports for it. 

So far what we have considered is the case that Quintessence is homogeneous in space. Since Quintessence is a 
scalar field one expects its fluctuations in space-time. In the following we study the isocurvature baryon number 

fluctuations generated in our model. _ 

Since the baryon number asymmetry is proportional to Q and in the tracking regime Q oc sjYiff), we have 


n_B 


oc 




(17) 


and 


o{ —) oc —, 

s 2VV 

S(nB/s) V' 

= w'*® ■ 

To get SQ, we are solving the equation of motion given below, 


(18) 


^Qk + ^HSQk H— K^Qk + V"5Qk — 0 , 

where SQk is the Fourier transform of SQ. To have an analytical formula for SQ we closely follow 
V as 


(19) 
by writing 


V — t^Hcq + t(1 ~ Cq){wb + Cq + 2)i7^ 


( 20 ) 
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where Cq is the sound speed of Quintessence defined by 


2 _ PQ 

= — = ’^Q - oTj/i I -T ■ 

PQ 3H{1 + wq) 

Consider the case that the sound speed is approximately constant [ 0 , the first term of the rhs. of Eq.(po|) can be 
dropped out, then the solution to Eq. ( 0 ) is plf : 

6Qk ^ r]~^f'^[CiJy{kri) + C2J-iy{kri)] , ( 22 ) 

where r] is the conformal time {dr] = ^), and a (x r] (x T~^ during the radiation dominated era. In Eq. (p^) , J±v{kri) 
are the first kind of Bessel functions, Ci, C 2 are constants which are time independent, and /3 are. 



-/3, 

~ Cq){u!b + Cq + 2) . 


(23) 


The modes of fluctuations which are interesting in cosmology are those outside the Hubble radius in very early 
universe, i.e. fcry <C 1. In such a limit the Bessel function can be approximated by 

Mkv) - 2^(1 + i.) ’ ^ 

On the initial value of SQk denoted by SQki, we assume the primordial fluctuation of the Quintessence scalar is 
generated during inflation. Being almost massless SQki is simply 

(25) 

where Hin is the Hubble constant during inflation. Furthermore the kinetic energy of the Quintessence will be diluted 
by inflation. Thus the Quintessence energy density remains potential energy dominated after inflation until the time 
rit when it goes into the tracking regime. This period of time is called the potential phase in the Ref. pl| . During this 
period pq is nearly constant and Cq ~ —wb — 2 ||^, hence /3 ~ 0 and V ~ 0. We can see that SQk is constant in 
the potential phase. But in the regime of tracking {rj > 774 ), 

Cq = wq , 

/3 = ^(1 - wq){wb +wq + 2) = ^(1 - wq)C- + wq) . (26) 

Numerically —1 < wq < 1/3, so we have (3 > 1/4. We obtain that 


<Q , 


V = i\iy\ 


J±u{k'n) 



exp(±j|i/| In ^) 

r(i±*kl) 


(27) 


Substituting J±„{kri) above and SQki (Eq.(|^) into (p^, we will get SQk- 
We have also numerically solved Eq. 0 ) for the Quintessence model we use in this paper. In Fig.l we plot SQk as 
a function of red-shift z. One can see that SQk is constant in the potential phase i.e. SQki = SQkt, then oscillates 
with amplitude which decays like oc 77 “^/^ during tracking. This confirms our analytical analysis above. 

Denoting the SQk at the temperature Tjj when the baryon number asymmetry freezes in by SQko, we have the 
ratio of \SQkD\ to the \SQki\ 


^1 <(-)“'’= 

SQki r]D It 


(28) 


where Tt is the temperature when Q transits from the potential phase to the tracking regime. Making use of the 
initial value of SQk generated during inflation, SQki = ^*(( 3 , we have 
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Following we calculate the RMS fluctuation of baryon number per logarithmic interval in our model which is 
given by 


{—)k<^\v\QD)IV{QD)r-;j^) 

riB 47r It 


Td-,1/2 


For the Quintessence model with potential in Eq.(] 


< 


A 




UB ' ^ 47r Mpi '■ Tt 

-Lt 


(30) 


(31) 


which will be smaller than 10 ® since Tt > Tq. To get the result above we have used the constraints on < 10 ® 
from the limit of the detecting of primordial gravitational waves predicted by inflation models. In addition we 
have also taken ~ 0(10) to satisfy the constraints from the nucleosynthesis. The isocurvature baryon number 
fluctuation in (|^ ) is small and is not conflict with the Boomerang and Maxima-1 data on the anisotropy of cosmic 
microwave background p^ . 

The isocurvature baryon number fluctuation in the spontaneous baryogenesis has been studied by Turner, Cohen 
and Kaplan where they And the amplitude sizable. In their model the potential of the ilion field 9 is rn?9^. The 
produced baryon asymmetry is linearly proportional to the initial value of the ilion held 9q. Therefore quantum 
fluctuations in the initial 9q during inflation will directly induce isocurvature baryon number fluctuations. In our 
model the Quintessence has the tracking behavior, and the produced baryon number asymmetry is almost independent 
of the initial conditions of Quintessence held. Consequently a small isocurvature perturbation is expectable. 

The tracking behavior of the Quintessence, however plays an important role in our discussions. Because of the 
tracking, Q can be large enough and the baryon number asymmetry required can be generated. For models without 
tracking property, such as the model with potential V{Q) = A^(l ± cos ^), we have checked and found that Q is too 
small to account for ub/s ~ 10“^°. 

Our mechanism for leptogenesis, however does put a lower limit on the reheating temperature Tr. It requires Tn 
be higher than Tt and Tq. For cut-off M of Eq.(l) in the range of 10^® ~ 10^®GeV, To is about 10^ ~ lO^^GeV, 
which will not put a strong constraint on models of inflation and reheating. 

In summary we have introduced in this paper an effective interaction of the Quintessence scalar to the ordinary 
matter and studied the possibilities of baryogenesis. Our model can naturally explain the baryon number asymmetry 
of our Universe ub/s ^ I0“^° via spontaneous leptogenesis. Our scenario provides a unified description for the present 
accelerated expansion and baryon number asymmetry generation of our Universe by the dynamics of Quintessence. 
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FIG. 1. Plot of SQk as function of red-shift 2 . The a;-axis is — ln (2 -|- 1), and the y-axis is \n\SQk/5Qki\- In this plot, Tt 
corresponds to temperature around 10^"* GeV and To — 10^° GeV. 





























